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Abstract 

Individual movements of a rod-like self-propelled particle on a flat substrate are quantified. 
Biological systems that fit into this description may be the Gram-negative delta-proteobacterium 
Myxococcus xanthus, Gram-negative bacterium Escherichia coli, and Mitochondria. There are also 
non-living analogues such as vibrated polar granulates and self-driven anisotropic colloidal particles. 
For that we study the Brownian motion of an asymmetric rod-like rigid particle self-propelled at 
a fixed speed along its long axis in two dimensions. The motion of such a particle in a uniform 
external potential field is also considered. The theoretical model presented here is anticipated to 
better describe individual cell motion as well as intracellular transport in 2D than previous models. 
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1. INTRODUCTION 



In classical physics the Brownian motion is the erratic motion of microscopic particles 
due to non-zero value of net stocastic forces at any given time exerted on them by atoms 
or moleculs of the surrounding medium. Such microscopic particles don't play active role in 
the motion and their irreversible dissipation of energy is compensated by reversible thermal 
fluctuations as expressed in the fluctuation-dissipation theorem. 

In contrast, if the particles have internal or external energy source leading to self- 
propulsion, their motion will have directional biasness. Such motion is reffered to as the 
active Brownian motion. Some example subjects of active Brownian motion relevant to our 
study are self-driven colloidal particles, living microorganisms, and vibrated Janus parti- 
cles [IH5]. These systems exhibit many interesting nonequilibrium effects understanding of 
which have been the research focus of many scientists in recent years. 

In this context, a phenomenological model of active Brownian motion has been devel- 
oped based on principle of energy take up, store and dissipation into kinetic degrees of 
freedom [HIE]. Other research team, J. R. Howse et ai, experimentally studied active Brow- 
nian motion of spherical polystyrene particles that are self-propelled due to asymmetric 
surface activity [8]. They interpreted the data using an expression for the mean-squared dis- 
placement that takes into account rotational diffusion of the self-propulsion velocity. More 
recently motion of deformable self-propelled particles analytically has been studied assum- 
ing that the particles can deform from circular shape when the propagation velocity is 
increased j9]. It was shown that such a particle undergo a bifurcation from straight motion 
to a circular motion by increasing the propagation velocity. Also, Brownian dynamics of a 
microswimmer composed of three spheres which propels itself in a viscous fluid by spinning 
motion of the spheres under zero net torque has been studied [10] . Finally, in ref [2] dynam- 
ics of a Brownian circle swimmer has been studied when the driving force does not coincide 
with the propagation direction. It was hypothesized that when the self-propulsion force and 
the particle orientation are in line, the motion would be along a straight line just perturbed 
by random fluctuations. 

In this paper, we present detailed study of the active Brownian motion of an asymmet- 
ric particle when the self-propulsion force and the particle orientation are in line. Unlike 
previous studies our approach takes into account the particle asymmetry. 
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FIG. 1: The diagram shows an ellipsoidal rigid particle self-propelled along its long axis. 

Collective behavior of such rods, which interact with each other through excluded volume, 
has been theoretically studied [11 J and it was shown that the self-propulsion enhances the 
longitudinal diffusion coefficient. Our calculations, however, show that in the case of one 
rod both the longitudinal and the transverse diffusion coefficients are enhanced and assume 
the same value at long times. 

The paper is organized as following: After introducing the model, we start by calculating 
the mean squared displacements and corresponding diffusion tensors for two cases: when 
a rod is self-propelled at constant speed and when such a self-propelled rod is in external 
constant potential field. The paper is concluded by the disucssion of the results. 

2. THE MODEL 

The model consists of an ellipsoidal rigid particle that is self-propelled along its long axis 
and undergoing Brownian motion in two dimensions (see FIG. py). Note that this model can 
be generalize to higher dimensions using the Pythagorean theorem. 

At a given time t the rod can be described by the position vector of its center of mass 
r(t) and the angle 9(t) between the x axis of the lab- frame and the long axis of the rod. In 
this frame the self-propulsion speed, which is taken along the long axis of the rod, is given 
by vn(t), where n(t) = (cos#(£), sin #(;£)) is a unit vector along the long axis of the rod. In 
the presence of external force F and torque r the dynamics of the rod is governed by the 
coupled Langevin equations, given by 




(1) 
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d t e = r R T + u(t), (2) 

where i — x,y for 2D and T^, T R are the mobilities relating velocity and angular velocity 
to force and torque, respectively. In the body-frame denoting the mobility along the long 
axis of the rod by T\\ and perpendicular by the mobility tensor in the lab-frame takes 
the following form 

Ar 

r ij (n)=F5 ij + —M ij (e), (3) 



cos 29 sin 29 



where Mij(9) = I I , T = F|l ^ r± , and Ar = T» — T±. The rotational noise 

sin 29 -cos 20 

£r is a Gaussian random variable with zero mean and variance 

(Ut)Ut')) = 2k B TT R 6(t-t'), (4) 

whereas £, x (t) and £ y (t) are Gaussian at a fixed 9(t) with zero mean and variances depending 
on orientation at a given time 

((MtjWUt) = 2k B TT ij (e(t))S(t - 1'), (5) 

where k R is the Boltzmann constant and T is the effective temperature. The ensemble 
average is calculated using (• ■ ■ ) = ^ "(' ' ' )do^o where (■ ■ • )g indicates average both 
over £j(£) and & f a fixed initial angle. 

The self-propulsion force doesn't induce torque, which in turn enables the rod to rotate 
freely about the axis passing through its center of mass and perpendicular to the plane 
of motion. This means that the angular displacements obey the Gaussian statistics with 
constant rotational diffusion coefficient. In the lab-frame the shape asymmetry of the rod 
makes the translational and rotational displacements coupled, which complicates the ana- 
lytic examination of the motion. The Gaussianity of the angular displacements, however, 
simplifies the matter and leads to detailed analytic results. 



3. SELF-PROPULSION AT CONSTANT SPEED 



As we mentioned previously, to carry out the calculations analytically we shall set r = 0. 
Therefore in the absence of external torque integration of Eq. ^ gives A9(t) = ^ R (t')dt' . 
Then with the aid of Eq. Q for the mean-squared angular displacement (the second moment) 
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we obtain 

(im)] 2 )= f [\ut , )Ut"))dt'dt" = 2D R t. 
Jo Jo 

Subsequantly the noise avarages of the sinusoidal functions will be evaluated using the 
following identity [T3] : 

^ e i[siA6»(t 1 )+s 2 A6»(t 2 )]\ )e _ e -D R [s 2 1 t 1 +slt2+2s 1 s 2 mm(t 1 ,t2)] ^ 

where D R = k R TT R is the rotational diffusion coefficent. 

Following the derivation of Ref. [H] in the absence of external force F one can integrate 
Eq. Q with respect to time and use Eq. (A6) to calculate the ensemble average of the 



sinusoidal functions. Then for the mean displacement one obtains 

(Ar} eo = UTi(cos0o,sin0 o ), (7) 

where r n = l ~ e n ^ Rt and t r = 1/2D R is the rotational time required for one radian diffusion. 
Finally, using Eq. for the magnitude of the mean displacement of the rod one has 

|(Ar) eo | = [{Ax)l + (AyJJJ* = 2vr R [l - e~M- (8) 

To find the mean-squared displacement tensor resulting from the integration of Eq. 
with F = 0, we note that on averaging over the noise, cross terms such as (6(^i) n j(^2)) 
vanish since they are first order in the noise. After doing so, in the expression for the 
mean-squared displacement tensor we are left with only two terms 

(Axi^Axjit))^ = I dh [ dt 2 (Uti)U t 2))e 
Jo Jo 

t ft 



+v 2 / dt x \ dt^it^n^))^. (9) 
Jo Jo 

The first term in Eq. (|9]) has been calculated in ref. [H] and the second term is given 



by Eq. (A7). If we put them together, for the mean-squared displacement tensor at a fixed 



initial orientation we obtain 



(AxiAx^ = [2Dt + 2r R v 2 (t - n)]5, 



'j 



+[ADt 4 + 2 ^-(r 1 - r 4 )}M t] (6 ), (10) 
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and the mean-squared displacement 



(Ar 2 ) = ADt + Ar R v 2 [t - 2r R (l - e"*/**)], (11) 
with D = k R TT and AD = k R TAT. For the displacement diffusion tensor of an active rod 



Eq. (10) yields 



D«{t, O ) = (AX f t Xj)e ° = [D S ~ TrV 2 ^ 



where Dg = D + t r v 2 . Averaging Eq. (12) over the initial angle, 6q, yields a time-dependent 
displacement diffusion tensor, 

D iM = 7T r ^oAi(*, ^o) = [D a ~ r R v 2 ^]5 ir (13) 

Z7T Jq t 

This shows that the system of non-interacting self-propelled rods is no longer diffusive 
as the displacement diffusion tensor is time dependent. However, in all cases the time 
dependence becomes negligible for t ^> tr, and as t tends to infinity unlike the passive 
ellipsoidal rigid particle case, for the self-propelled ellipsoid the single long-time diffusion 
coefficient is enhanced according to 

-> DsSij. (14) 

When t 3> t r the directional memory is lost and the rod behaves as a perfect spherical 
Brownian particle with enhanced diffusion coefficient. 

On the other hand, when t <C t r the diffusion tensor and the mean-squared displacement 
are give by D^t) ~ [D + ^]<% and (Ar 2 ) ~ ADt + v 2 t 2 , respectively, indicating ballistic 
behavior at short times. 



4. SELF-PROPULSION IN EXTERNAL CONSTANT FORCE POTENTIAL 
FIELD 

In the previous section we studied the active Brownian motion of a rigid rod in the 
absence of external field and torque. Here we discuss the motion under the infulance of 
external constant force. As in the previous case we assume the rod is constrained to move 
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in a plane. An example of such a system may be active colloidal suspension of asymmetric 
particles under gravity. 

Without loosing any generality, we can always rotate the coordinate system (lab-frame) 
such that the external force is along the x axis, F = Fx. Then using Eq. @ from Eq. Q 
we can write 

1 r* 

Ax(t) = FTt + -FAT / cos 26dt' 

2 Jo 

+ / ixdt' + v I cos 9dt', (15) 
Jo Jo 

Ay(t) = -FAT [ sm20dt' + [ £ y dt' + v [ sm9dt', (16) 
2 Jo Jo Jo 

with averages given by 

(Ax) 9o = FTt + vti cos O + ^FAIY 4 cos 29 , 

(Ay) eo = vti sin 6 + -FAIY 4 sin 29 . (17) 

The averages of the mean-squared displacements along the x and y axes are given by 

(Ax 2 ) d0 = ^F 2 Ar 2 [3(t - r 4 ) + (r 4 - r 16 ) cos4# ] 
+ ADt a cos 26 + F 2 TATtr A cos 28 + 2vFTtr 1 cos 6 
+2Dt + F 2 T 2 t 2 + ^T R v 2 [3{t - n) + (ri - r 4 ) cos2# ] 

'" r^FAr[40(4r 1 - 3te~ Dfl * - r 4 ) cos# 



120 

+3(5ti + 8r 4 — 13r 9 ) cos3#o 



(18) 



(Ay 2 ) 9o = ^F 2 AT 2 [3(t - r 4 ) - (r 4 - r 16 ) cos40 o ] 
2 

-ADr 4 cos29 + -r R v 2 [3(t - n) - (n - r 4 ) cos26» ] 

+— r^FAr[40(4n - 3te-° Rt - r 4 ) cos# 
-3(5n + 8r 4 - 13r 9 ) cos3# ] + 2L>t, 



(19) 



7 



and the average of the cross-correlation between displacements along the x and y axes is 

(AxAy)e = ^F 2 Ar 2 (r 4 - r 16 ) sin4# + vFTtn x 

1 2 
sin 6^0 + -F 2 rArtr 4 sin 26»o + -trv 2 ^ - r 4 ) sin20 o 

+^r^FAr(5n + 8r 4 - 13r 9 ) sin30 o , 

(20) 

In our calculations we took into account that the terms odd in the noise average to zero. 
Then the mean-squared displacement is 

2X _ 1 (** M Moo + (^k 



i r 



2t 

= ADt + Ar R v 2 (t - n) + F 2 r 2 t 2 

+ ^F 2 AT 2 r R (t-n), (21) 

and the long time t 3> tr diffusion tensor averaged over the initial angle, 9q, yields the 
following experssion 

D«(t) = ± f* d6 {Ax ^ )so = D SF 5 l3 (22) 



with 



D SF = D + t r v 2 + vFATt r 

+—F 2 AT 2 t r + -F 2 T 2 t. (23) 
16 4 K J 

Self- propulsion vs. external potential field - We saw that for the self-propelled 
particle the long time single diffusion coefficient is enhanced according to D$ — D + r R v 2 . 
This enhancement is similar to that of a non self-propelled ellipsoid due to external constant 
potential field (Dp = D + r^|F| 2 Ar 2 /16 where |F| is the force magnitude of the external 
potential field [TH] ) since both of them are proportional to the square of the force magnitude. 
However, the latter explicitly depends on the asymmetry of the particle, Ar, but the former 
does not. The origin of this discrepancy may be apparent from the fact that our model 
assumes that in addition to the external constant force the rod has a fixed speed v directed 
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along its long axis at all times while the former model assumes the rod feels only an external 
constant force F at all times. Therefore similar to the hydrodynamic drag force, the external 
potential field brings in the friction factor between the two situations. 

Also, a clear difference between the two situations is easily understood by neglecting the 
noise and hence considering a deterministic description. Then one can immediately calculate 
the velocity and displacement in the lab frame as functions of time and the initial angle, 
6q, for given initial conditions. In the absence of external potential field we have zero mean 
displacement and zero velocity when averaging over 8q. However, in the presence of external 
potential field one obtains a mean velocity TF X in the x direction and TF y in the y direction 
with corresponding mean displacements. 

5. DISCUSSION 

In writing down the Eqs. ([TJ and ^ we have dropped the inertial terms assuming the 
motion is overdamped, that is, the viscous terms dominate the inertial ones. If we had kept 
those terms, they would have generated homogenous part of the solutions which exponen- 
tially decay with the characteristic relaxation times, m x for translational and I x Tr 
for rotational motion. Here m is the mass and / is the moment of inertia about the z axis 
passing through the center of mass of the rod. Our results are, therefore, not valid in time 
intervals comparable to the relaxation times given the fact that recently developed exper- 
imentation technics, namely, use of the Photonic Force Microscopy allows high precision 
measurements in such small time intervals [T2] . 

Mean displacement - Eq. ^ has limiting forms of |(Ar)# | = 2vtr for t ^> tr sug- 
gesting that at very long times the self-propelled rod must behave as a spherical Brownian 
particle in the region centered Ivtr far from its initial position and |(Ar)e | = vt for t <C tr. 
The latter is intuitive since at early times the motion is anisotropic which is not the case 
for the long time result. To give a full physical meaning to these results, we contrast the 
cases when the self-propulsion (SP) is 'turned on' and 'off'. First thing to note about these 
two situations is that both of them have the same rotational time, tr, since the net torque 
is zero in both cases. Using this, we can picture what the rod would do before and after 
the crossover has occurred when the SP is on and when it is off: (i) before the crossover 
(t < tr), short time anisotropy gets contribution from SP which causes the rod move faster 
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in a particluar direction than when SP is turned off. Therefore in the anisotropic regime the 
rod displaces more due to SP than when SP is turned off, (ii) after the crossover has occurred 
(t > Tr)> m both cases the rod starts performing isotropic diffusion about the position that 
it has arrived during the time while t < tr. Therefore 2vtr is the mean distance between 
the centers of isotropic diffusion domains of non self-propelled and self-propelled rods. 

Note that if we average over the initial orientation, components of the mean displacement 
are zero and so is the mean displacement. This is not in contradiction with what we have 
above. Averaging over the initial angle is equivalent to averaging over the number of non- 
interacting self-propelled rods that started up from the same point in isotropic manner. 
At long times, the centers of isotropic diffusion for these rods are distributed isotropically 
about the starting point and hence have zero mean displacement. From this point of view, 
we average the diffusion tensor over the initial orientation to analyze transient behavior of 
diffusion in the system of non-interacting such rods. 

A particle undergoing classical Brownian motion takes the same amount of time for going 
from a given initial point to a final point and coming back. For our model this is not the 
case: The domain of isotropic motion is centered 2v tr far from an arbitrary starting point 
implying that a particle would need more time (equivalently less probability) for coming 
back. In fact, it has been observed that when the pigment granules are spreading out from 
nucleus of the cell, they go straight through actin filament intersections until they get to 
the end of the filament, but if they need to return to the nucleus, cargos often switch at 
actin filament intersections hence taking relatively longer time to come back [T7] . This is in 
accordance with the prediction of the model presented here. 



Mean-squared displacement - In Eq. (11) the second term is exactly the same as 
the one in refs [18-20J obtained for the persistent random walk [21] with persistence time 
P = 2tr and the weak noise limit of the expression obtained from a model of Brownian 
particles which are pumped with energy by means of a non-linear friction function (upon 
setting tr = Vq/4D) in refs [6] [7], while the first term is the contribution due to classical 
Brownian motion. This is similar to that of the Levy-walk superdiffusive motion which can 
be decomposed into thermal orientation fluctuations and an active motion of the rods with 
a constant velocity along their long axis [31]. 

It is known that at short times the cell appears to persist at the direction of motion but at 
long times the pattern of movement is similar to the classical Brownian motion [22] • Because 
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of this, the persistent random walk model have been used to quantify cell migrations [23- 
The results presented here suggest that by including contribution due to the classical 



Brownian motion as in Eq. (11) one can describe individual cell movements on a plane 
substrate more accurately than it has been using the persistent random walk model. 

Recently, an experiment has been carried out by Jonathan R. Howse, et al. who studied 
the motion of polystyrene spheres with a diameter of 1.62 /im immersed in hydrogen peroxide 
and half-coated with platinum [Sj. Due to the local osmotic pressure gradient created by 
the asymmetric chemical reaction such spheres were propelled. They were able to fit the 
experimental data only with the short- and long-time limits of the expression for their model, 



which is the same as Eq. (11 ) if we replace tr with Tr/4. But if we consider the displacement 
components along the x and y axes separately, then immediately we see a major difference 
between these two cases which has been hidden in the expressions for the mean-squared 
displacements. The first model consists of a self-propelled sphere whose direction of motion 
is subject to diffusion and hence no expression could contain shape asymmetry as it is zero. 
In contrast, our model contains shape asymmetry, AD, which makes it distinct from any 
previously studied model when considering the x and y components of quantities of our 
interest, such as the mean-squared displacement. We are able to recover all previously 
obtained results upon setting AD = in our calculations. 

In summery, here we present a simple model to describe the individual movements of 
rod-like living microorganisms on a plain. For that we studied the Brownian motion of 
an asymmetric rigid particle self-propelled at a fixed speed v along its long axis in two 
dimensions. The results for the case when such a particle is in external potential field are 
also presented. Some relevant biological systems are rod like living microorganisms such 
as the Gram-negative delta-proteobacterium Myxococcus xanthus, Gram-negative bacterium 
Escherichia coli and Mitochondria. Non-living self-propelled systems which resemble our 
model may be a polar rod on a vibrating substrate pQ, self-driven anisotropic colloidal 
particles and catalytically driven nanorods [21 E]. These results are expected to better 
quantify motion of such systems over the entire range of time. 
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Appendix A 

Here all averages (...) are taken at a fixed initial angle, 8 , and 9\ = #(£2) and 0i = #(£2) 
We start by writing 



(7li(*i)7lj(* 2 )) 



(cos 9i cos 62) (cos 6*i sin 2 ) 
(sin 6*i cos 62) (sin 6\ sin 2 ) 



(Al) 



and 



cos Q\ cos 6 2 = - [cos(#i — 6 2 ) + cos(0i + 6 2 )] j 
cos 0! sin 9 2 = ^[sin(#i + 2 ) - sin(0i - 2 )], 
sin Q x cos 2 = 2 t sin (^i - ^2) + sin(0i + 2 )] , 
sin0! sin 2 = ^[ cos (^i — $2) _ cos(0! + 2 )]. 



(A2) 



To calculate the averages we use Eq. (|9p and the fact that cos0 = $te l9 and sin0 = ^e 19 
where 9ft and 3 stand for the real and imaginary parts respectively. Then 

2(cos0(ti)cos0(t 2 )) = e - DR[tl+t *- 2m[n{tl ' t2)] 
+ cos29 e- DR[h+t2+2m ' m(tl ' t2)] , 

2(sin0(ti)sin0(t 2 )) = e -DR[ti+t2-2mm( tl ,t 2 )] 

-cos20 o e- D «[* 1+ * 2+2min ^'* 2 ^ 
2(cos0(t 1 )sin0(t 2 )> = s in20 o e- Dii[ * 1+t2+2min( * 1 '* 2)] , 
2(sin0(t 1 )cos0(t 2 )) = sin20 o e- Dfl[ * 1+t2+2min( * 1 '* 2)1 . 

(A3) 



12 



If we put Eq. (A3) into Eq. (Al), we obtain 



+M ij (6 )e- DR[tl+t2+2min{tut2)] 



Now we need to integrate Eq. (A4) twice. To do that we use the following formulae 



t pt pt rt\ pt rt2 

dti / dt 2 ... = / dti / dt 2 ...+ / dt 2 / dti..., 



o jo 



o jo 



ti > t 2 



o Jo 



ti < t 2 

which can easily be proved using geometrical considerations. Then we have 



f dh f dt2e- D ^ +t * +2min ^^ 
Jo Jo 

2 f dti [ tl dt 2 e- D ^ +3 ^ = ^[n(t) -n(t)), 

Jo Jo 

f dti r^ 2e -^[*i+^-2min( tll t 2 )] 

Jo Jo 

= 2 f dh / 1 dt 2 e- DR[tl - t2] = 4r R [t - n(t)]. 
Jo Jo 



Finally, if we put Eqs. (A4), and (A6) together, we obtain 



t rt 



v dti dt 2 {ni(t 1 )n j (t 2 ))e 
Jo Jo 
2 

2r R v 2 [t - Ti}5ij + -t r v 2 [ti - r 4 ]My(0o)- 



(A4) 



(A5) 



(A6) 



(A7) 
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